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ABSTRACT 

The interaction between the tire and road has long been of 
interest for vehicle dynamic simulation. A planar tire model is 
developed to capture the tire circumferential displacements and 
calculate the spindle force according to the tire shape. The tire 
is discretized into segments and Hamilton s principle is used to 
derive the model mathematical expression. It is shown that the 
static constraint modes are functions of two non-dimensional 
parameters; a third parameter defines the overall stiffness. 
These parameters are experimentally identified for a specific 
tire. The bridging and enveloping properties are examined 
circumferentially. The prediction accuracy of spindle force 
with respect to tire-road interference is evaluated by comparing 
the simulation and experimental response for a quasi-static 
cleat test. The simulation result of spindle force agrees with the 
experimental data and the process can be implemented as a 
morphological pre-filter of road profiles for more efficient 
vehicle modeling and simulation. 
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NOMENCLATURE 

k Distributed spring rate. 

u Radial deflection. 

R Tire radius. 

x F ac Vector of active constraint degrees of freedom. 
O c Matrix of constraint modes. 

F A q Force acting on one single tire piece, 
f Linear distributed load density. 

1 length of one single tire piece. 

A 6 Angular range of each tire segment. 

V b Bending potential energy for each tire segment. 
o xx Elastic stress for each tire segment. 

£ xx Elastic strain for each tire segment. 
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M Bending moment for each tire segment. 

E Elastic modulus. 

/ Second moment of area for each tire segment, 
p Massive density. 

Y Shear strain. 

V s Shear potential energy for each tire segment. 
G Shear modulus. 

T Kinetic energy for each tire segment. 

V e Elastic potential energy for each tire segment. 
W External work on each tire segment, 
n Tire segment number. 


INTRODUCTION 

The interaction between the tire and road has long been of 
interest for vehicle dynamic simulation and chassis design. 
Therefore, a proper tire model that captures the circumferential 
deformation and predicts the load is required. The tire model 
should bridge cracks in the pavement and envelope small 
bumps, and the resultant force should be obtainable through the 
circumferential deformation. A considerable amount of 
research has been done in this area. Umsrithong presents a 
detailed tire model review in his dissertation!!]. Among all 
kinds of tire models, the in-plane ring model with elastic 
foundation places most emphasis on the circumferential 
deformation. Badalamenti et al. developed the radial spring tire 
model, which makes the radial spring element deflection 
dependent on the adjacent element deflections[2]. This model is 
shown to accurately describe the tire envelopment behavior. 
Zegelaar and Pacejka used a flexible ring model to simulate the 
quasi-static response of a tire rolling over an uneven surface[3]. 
The simulation results agreed with experimental data, but the 
complexity of the model increases the computational 
requirements. Loo modeled the tire by using a flexible circular 
ring under tension with a nest of linear springs and dampers 
arranged radially[4]. The ring tension and foundation stiffness 
are obtained by experimental data from a static point-load test 
and the contact patch measurements. This model takes tire 
rolling resistance into consideration. Gillespie used a simple 
radial spring model to simulate the stiffness variation 
circumferentially [5]. The magnitude of the radial force 
variation is found to be relatively independent of speed. 
Takayama et al. developed a mass-spring tire model to predict 
the transient response of a tire encountering a cleat [6]. The 
belt and tread region was modeled by a rigid ring and 
deflections from the cleat are absorbed by a linear and planar 
spring attached to the rigid ring. 


Ferris proposed the static constraint mode tire model to 
capture the tire enveloping properties [7]. The model applies an 
iterative process to identify the active constraints at the tire- 
road interface from which the resulting circumferential 
deformation of the tire is formulated. This work developed in 
this paper improves on this static constraint mode model. The 
theory is expanded to include shear effects and the required 
stiffness parameters are reduced from five dependent 
parameters to two parameters that define the constraint modes 
and an independent parameter that defines the overall tire 
stiffness. In this way the parameters that define the shape of 
the tire are independent from the overall stiffness. The ability 
of the model to generate the correct tire shape is validated for 
cracks and bumps to check the enveloping and bridging 
properties. The appropriate overall tire stiffness parameter is 
identified through the experimental result from a single flat- 
plate test. The resultant spindle force can be calculated through 
the circumferential displacement. A set of experimental quasi¬ 
static cleat tests are used to validate the entire modeling 
process. 

The work is organized as follows: First the background on 
the constraint mode method is reviewed. Next the improved 
theoretical deduction process of the constraint mode model is 
developed. The model parameters are then identified via the 
deformed tire shape and validated with experimental data. The 
bridging and enveloping properties of the model are examined 
to check the effectiveness of the model in terms of 
circumferential deformation. Finally, the simulation result of 
spindle force versus deformation over a cleat is compared with 
experimental data to demonstrate the force prediction accuracy. 


BACKGROUND 

Since the development of component mode synthesis 
originated by Hurty and Gladwell, the reduced representation of 
component model has long been of interest in developing 
component constraint modes [8, 9]. Guyan provides a static 
reduction method for mass and stiffness matrices [10], from 
which Kuhar and Stahle developed the dynamic transformation 
method to derive constraint modes from the response at a given 
frequency [11]. Ferris proposed a constraint mode tire model to 
describe the static deflection where the system is assumed to be 
adequately represented by an undamped, second order, linear, 
time invariant, ordinary differential equation as given in Eq. 
( 1 ). 

[m]{u} + [k]{u} = {f} (1) 

The physical degrees of freedom can be categorized into 
active and omitted ones, represented by ‘a’ and ‘o’ respectively. 
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In general, the boundary degrees of freedom must be a subset 
of the active degrees of freedom. In the tire modeling case, the 
active degrees of freedom must include those in contact with 
the terrain. Therefore the mass and stiffness matrices are 
reordered and partitioned as given in Eq. (2). 


w -[£2 !' 


[m aa ] 

[m ao ] 

[m oa ] 

[m 00 ] 

■[ k aa ] 

[k ao ] 

[k oa ] 

[k 00 ]. 


( 2 ) 


By substituting Equation 2 into Equation 1, and following 
the work by Kuhar and Stahle, the constraint frequency Gl>c and 
the constraint modes can be computed as shown in Eq. (3). 


(\[k aa ] 

vL[fc oa ] 



[m ao ]l\ r{H» ac }'l 

[m 00 ]JJl{V 0C }j 


(3) 


The lower set of equations can be used to derive the 
constraint mode matrix given in Eq. (4). 



FIGURE 1. IN-PLANE TIRE MODEL SCHEME 



({’P ac }| 

( {V oc } j 

= [ [/] 

L-([fc 00 ] - a>c [m 00 ]) - 1 ([fc 0a ] - a) c 2 [tn oa ]) 

= [0 c ]*F ac 


{ipac} 


(4) 


The active constraint degrees of freedom { l F ac } are defined 
by iteratively applying the unilateral geometric constraints at 
the tire-road interface [7]. 

The tire is assumed to be axisymmetric and isotropic 
circumferentially. A simple model which captures the essential 
in-plane tire behavior is shown as Figure 1. The tire is modeled 
as an inextensible ring supported by an elastic foundation 
representing the radial stiffness. The tire is divided into n 
segments and analyzed discretely as Figure 2 shows. Each 
individual tire segment is considered as an Euler elastic beam, 
{u} is the radial deflection used to describe the circumferential 
tire displacement. A recursive algorithm to determine the 
active geometric constraint is developed in Ferris’ work [7]. 


FIGURE 2. TIRE RING SEGMENTS IN 3D 


MODEL DEVELOPMENT 

In order to examine the relationship between the stiffness 
matrix and the tire properties, a more detailed expression of the 
model is required. The shear potential energy term is added to 
Ferris’s constraint mode deduction [7] to describe the relative 
movement of two adjacent tire segments. Therefore the 
equation resulting from the application of Hamilton’s principle 
for one tire segment can be rewritten as Eq. (5). 

f t2 (5) 

(6T - SV b - SV e - SV s + 5W)dt 

K 

= 0 Vt 1# t 2 

The kinetic energy term vanishes under the quasi-static 
assumption. Each term and the corresponding virtual changes 
are computed as Eq. (6) - Eq. (10) : 
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1 (du\ 2 

T= z pi y =° 


V, 


Vs = 


_ Ell /3 2 u\ 

2 \dx 2 / 

1 

V = — klu 2 
e 2 

Gblh /du\ 2 


2 \dx) 


W = flu 


5T = 0 

(6) 

2kTiR 

C = kl = 

n 

Elastic stiffness parameter 

3 4 u 

§V b = EIl —6u 

(7) 

2fnR 

f a 0 - n 

Force on each tire segment 

SV e = kludu 

(8) 



d 2 u 

5V S = Glbh—r6u 
dx 2 

5W = fl6u 

(9) 

(10) 

In this work, the system is assumed to be quasi-static, therefore 
the constraint frequency is zero and Eq. (4) reduces to the 
Guyan form where the circumferential displacements of the tire 
can be predicted by Eq. (16). 


The equation of motion for each segment is then simplified to 
Eq. (11). 



[I] 

— ([k 00 ]) -1 [k oa ] 


j"vp ac } _ [0 c ]{ x F ac } 


(16) 


d 4 u d 2 u (11) 

Ell —— + klu + Glbh ——- = fl 
dx 4 dx 2 


The solution to Eq. (11), assuming that the solution is 
separable in space and time, can be expressed as Eq. (12). 

U( x ,t) = V x e At (12) 


The stiffness matrix [k] has the following form. 


1 

(Xi 

a 2 

0 • 

• a 2 

(xt 


Oti 

1 

(Xi 

ot 2 • 

• 0 

(X 2 


ot 2 

0 

(Xi 

ot 2 

1 

(Xi 

(Xi 

1 • 

0 

•• 0 

0 

0 

= k 0 [a] 

_0Ci 

a 2 

0 

0 • 

•• 

1 . 



(17) 


Therefore the fourth derivative with respected to x can be 
expressed as Eq. (13). 



At 


(13) 


The derivatives of V 0 with respect to x at point zero, i.e., 
the derivative of u with respect to x, can be approximated by 
adjacent points with finite difference method. Eq. (14) shows 
the approximation of V 0 with 2 adjacent points on each side. 

„„ U_ 2 - 4u_! + 6u 0 - 4u x + u 2 (14) 

U 0 - 14 


By substituting Eq.(14), Eq. (11) is transformed into Eq. 
(15). 

Au n _! + (B - 4A)u n + (6A - 2B + C)u 0 + (B (15) 
- 4A)u x + Au 2 = F A0 


Where k 0 = 6A — 2B + C, a x = B — 4A, a 2 = A 

Eq. (17) can be reordered as: 

r a ao a aa j r u a | _ rF aa | 

0 La oa a 00 -! lu°l l 0 J 


(18) 


It can be seen that the resultant tire shape depends 
on a 00 and a oa , therefore the parameters a x and a 2 completely 
determine the constraint modes and together with the active 
degrees of freedom completely specify the radial tire 
deformation. The parameter k 0 is an overall stiffness 
coefficient used to tune a specific tire. The constraint model 
parameters must be estimated and adjusted according to 
mechanical properties of different tires. A parameter 
identification process is necessary so that the model can be 
applied to force predictions of different tires with a consistent 
practice of identifying the set of modal parameter values. 


Where 

El Eln 3 Bending stiffness parameter 

¥ = (2ttR) 3 

Gbh Gbhn Shear stiffness parameter 

” ~T~ 2ttR 


PARAMETER IDENTIFICATION PROCESS 

The experimental data used throughout this work is 
provided by Professor Schalk Els of the University of Pretoria, 
South Africa. The tire used is Continental Conti-Trac AT 
238/85 R16 and the schematic experiment set up is shown as 
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Figure 4. During the Parameter Identification Process, a fiat 
plate is pushed towards the center of the tire gradually and the 
resultant spindle force F is measured with respect to the tire- 
plate interference e. No data was available for the deformed 
shape of the tire, so the parameters a x and a 2 are estimated by 
trial and error practice so that the tire maintains a reasonable 
deformed shape. The results of this choice of shape parameters 
are developed in the Validation section of this work. 



An optimization process with respect to the difference 
between the simulated and the experimental forces is employed 
to identify the overall stiffness parameter, k 0 , for a specific tire. 
The comparison between the simulated result and the 
experimental data shows very good agreement as indicated in 
Figure 5. The model parameters for the specific tire used in 
this work given the number of finite segments, n, is provided in 
Table 1. There were 72 segments used in this work. 



FIGURE 5. COMPARISON BETWEEN THE FLAT-PLATE 
EXPERIMENTAL AND SIMULATED FORCE, F, WITH 
RESPECT TO THE INTERFERENCE, E. 


TABLE 1. TIRE SHAPE AND STIFFNESS PARAMETERS FOR 
A PARTICULAR TIRE 


Parameters 

Values 

a ± 

1388.8889n — 12.8601n 3 

a 2 

3.2150n 3 

k 0 

19.2301n 3 - 2777.7778n 

1 

+ 5400000- 
n 


VALIDATION 

The tire model can be tested over small cracks and bumps 
using the tire shape parameters, a x and a 2 , listed in Table 1. 
The resultant shapes are shown in Figure 6. The black and the 
red lines represent the original and deformed tire shape 
respectively, while the blue lines represent the discretized 
terrain profile. The tire bridging and enveloping properties are 
evident. 
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FIGURE 6. BRIDGING AND ENVELOPING PROPERTIES OF 
CONSTRAINT MODE TIRE MODEL 

A cleat test is conducted using the test layout as shown in 
Figure 7. A flat plate with a square tube at the center is pushed 
gradually towards of the tire and the resultant spindle force F is 
measured with respect to the tire-plate interference e. The side 
length of the square tube is 19 mm. 



The comparison between the simulation and measured data 
shows reasonable agreement in linear region, i.e. the small 
interference region (up to 0.05m), as shown in Figure 8. It 
appears that the stiffness parameter, k 0 , that was identified 
using a flat plate underestimates the stiffness for a cleat test. 
UNCLASSIFIED: Dist A. Approved for public release 


This may be due to the fact that the shape of the tire during the 
flat plate test is unknown. A better estimate of the tire shape 
parameters, a x and a 2 , may in turn lead to a better estimate of 
the stiffness parameter. The force versus deflection curve as 
shown in Figure 8 shows a discontinuity in slope when the 
deflection is around 0.047 m. In the first segment with smaller 
slope, the tire is suspended by the cleat event in the air and the 
contact area related to force generation is constrained to the 
cleat surface only. The ground does not provide any support at 
this point. Therefore the force change rate with respect to 
deflection is relatively small. In the second segment with a 
larger slope, the tire tread outside the cleat contact area touches 
the ground and the contact area increases. The ground starts to 
provide vertical support force to the tire. Therefore the rate at 
which the generated force changes with respect to the 
deflection increases. 



FIGURE 8. COMPARISON BETWEEN THE CLEAT 
EXPERIMENTAL AND SIMULATED FORCE, F, WITH 
RESPECT TO THE INTERFERENCE, E. 


DISCUSSION 

The constraint mode model can serve as a morphological 
pre-filter of actual terrain since it incorporates the enveloping 
and bridging property. That is, the wheel trajectory can be 
computed given a section of actual terrain profile over which 
the wheel travels and a constant spindle load. Figure 9 shows 
the comparison between the actual terrain profile and the 
adjusted wheel center trajectory. The wheel center height 
trajectory is vertically adjusted to compare with the original 
terrain profile. The terrain profile is smoothed due to the tire 
bridging and enveloping properties with respect to the 
undulations. For example, Figure 9 demonstrates that deep 
narrow cracks are filtered out of the resulting profile. In future 
work, a quantitated parameter identification process will be 

































developed to determine the deformed shape parameters a x and 
a 2 quantitatively instead of trial and error. The tire constraint 
model will be further improved so that it can be applied on not 
only deterministic road profiles, but also to stochastic terrain 
profiles and ultimately to deformable terrain. 



FIGURE 9. CONSTRAINT MODE FILTER. 


CONCLUSIONS 

The main contributions of this work are the enhancement 
to the constraint mode formulation for the tire-terrain interface 
and the development of the two tire shape parameters, 
a x and a 2 , and an independent tire stiffness parameter, k 0 . 
With this detailed expression of the model parameters, a model 
parameter identification process for tire constraint mode is 
established using only the flat plate results. A table with typical 
tire shape and stiffness parameters is provided. The model is 
shown to have bridging and enveloping properties over cracks 
and bumps, respectively. Reasonable predictions for the 
spindle force during a cleat test are demonstrated and the 
discontinuity in the stiffness during the cleat test is properly 
identified. This planar model can be used to predict the tire 
shape and the force acting on the spindle given a specific tire 
and a section of terrain profile. It can also serve as a 
morphological pre-filter of actual terrain to obtain the trajectory 
of the wheel center given constant load acting on the spindle. 
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